Packing Confined Hard Spheres Denser with Adaptive Prism Phases 
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We show that hard spheres confined between two parallel hard plates pack denser with periodic 
adaptive prismatic structures which are composed of alternating prisms of spheres. The internal 
structure of the prisms adapts to the slit height which results in close packings for a range of plate 
separations, just above the distance where three intersecting square layers fit exactly between the 
plates. The adaptive prism phases are also observed in real-space experiments on confined sterically 
stabilized colloids and in Monte Carlo simulations at finite pressure. 
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How to pack the largest number of hard objects in a 
given volume is a classic optimization problem in pure 
geometry [l[ . The close- packed structures obtained from 
such optimizations are also pivotal in understanding the 
basic physical mechanisms behind freezing 0, Q and 
glass formation Moreover, close-packed structures are 
highly relevant to numerous applications ranging from 
packaging macroscopic bodies and granulates |5| to the 
self-assembly of colloidal and biological 0,0] soft mat- 
ter. For the case of hard spheres, Kepler conjectured 
that the highest-packing density should be that of a peri- 
odic face-centered-cubic (fee) lattice composed of stacked 
hexagonal layers; it took until 2005 for a strict mathe- 
matical proof More recent studies on close packing 
concern either non-spherical hard objects (Io| such as el- 
lipsoids fill [lp . convex polyhedra [H Q ( in particular 
tetrahedra |l5j). and irregular non-convex bodies (l6j or 
hard spheres confined in hard containers (l7l - fl9| or other 
complex environments. 

If hard spheres of diameter a are confined between two 
hard parallel plates of distance H, as schematically illus- 
trated in Fig. [U the close-packed volume fraction </> and 
its associated structure depend on the ratio H/a. Typ- 
ically, the complexity of the observed phases increases 
tremendously on confining the system. Parallel slices 
from the fee bulk crystal are only close-packed for cer- 
tain values of H/a: A stack of n hexagonal (square) lay- 
ers aligned with the walls, denoted by n,A (nO), is best- 
packed at the plate separation H n & (-ff n rj) where the lay- 
ers exactly fit between the walls. Clearly, for the minimal 
plate distance H = ifiA = packing by a hexagonal 
monolayer is optimal. Increasing H/a up to H2A, a buck- 
led monolayer [2(J and then a rhombic bilayer 21 1 be- 
come close-packed. However, for H2A < H < H4A, the 



H 




FIG. 1. Schematic illustration of hard spheres of diameter a 
confined between two parallel hard plates of separation H. 



close-packed structures are much more complex and still 
debated. Both, prism phases with alternating parallel 
prism-like arrays composed of hexagonal and square base 



23j and morphologies derived from the hexagonal 



25| were proposed as 



close-packed (hep) structure [2J 
possible candidates. 

For confined hard spheres, the knowledge and control 
over the close-packed configuration is of central relevance 
for at least two reasons: First, the hard sphere system 
away from close-packing is of fundamental interest as a 
quasi-two-dimensional statistical mechanics model. At 
low densities, a hard sphere gas is stable, which will crys- 
tallize as the density is increased beyond some threshold 
value. As such, the model represents a classical route 
to understand freezing between two and three spatial di- 
mensions (26|. The associated fluid-solid transition will 
be strongly affected by the close-packed structure. Sec- 
ond, the confined hard sphere model is almost perfectly 
realized in nature by mesoscopic sterically-stabilized col- 
loidal suspensions [22j,|27| which can be confined between 
glass plates providing a slit-like confinement. At high im- 
posed pressures, colloids will self-assemble into the close- 
packed structures. It has been shown that this is the 
key for the controlled fabrication of nano-sieves and of 
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membranes with desired morphology [281 ] . 

In this Letter, we explore the close-packed structures of 
confined hard spheres by combining numerical optimiza- 
tion, experiments and computer simulation. Using a sys- 
tematic penalty optimization method, we find the whole 
cascade of close-packed structures in the range of plate 
distances -Hi a < H < i?4A- As an important building 
block for close-packing, an adaptive prism is identified 
which adjusts its internal structure flexibly to the slit 
height H/a. This prism has a base with a rhombic sym- 
metry and neighboring prismatic arrays are shifted rela- 
tive to each other. The resulting adaptive structure max- 
imizes the packing fraction in the regime beyond H 3 q. 
We also propose a further close-packing prism phase of 
square symmetry that packs densest in the regime just 
beyond H40 and shows a two-dimensional relative lateral 
shift between the prisms. We confirm the stability of the 
new adaptive prismatic structures both in real-space ex- 
periments on confined sterically stabilized colloids and in 
Monte Carlo simulations at finite pressure. In the follow- 
ing, we first describe the results from the penalty method, 
then discuss real-space data for confined colloidal samples 
and subsequently turn to Monte Carlo simulation results. 
Details of the numerical methods, experiments and sim- 
ulations are listed in the Supplemental Material [2!| . 

In our numerical calculations, we consider periodic 
structures with up to f2 particles per unit cell thereby 
covering all hitherto proposed structures j29[. To maxi- 
mize the packing fraction 0, we optimized the cell shape 
and the particle coordinates of these structures. How- 
ever, investigating the dense-packing of hard spheres ac- 
commodates a constrained optimization: the free vol- 
ume must be minimized under the constraint of non- 
overlapping spheres. To circumvent the discontinu- 
ous, constrained optimization, we employed the penalty 
method [3Cj | in our numerical calculations. By adding 
a penalty term in case the spheres intersect which de- 
pends continuously on the overlap volume, we obtained 
a continuous and unconstrained penalty function which 
can be minimized in the classic way to predict the opti- 
mal particle coordinates. The penalty method offers the 
flexibility to use a relatively broad range of candidate 
crystalline lattices and has recently been shown to allow 
a very efficient handling of packing problems [3l| . 

The resulting volume fractions of the densest packed 
phases are shown in Fig. [2] as a function of H/a in the 
regime between the hexagonal monolayer I A (H/a = I) 
and the triangular tetralayer 4A (H/a = V6 + I). For 
H 1A < H < H 2A the classic sequence [H ED, [IH f A 
£?— >2D— >2i?— >2Ais confirmed. Here, B is a buckled 
hexagonal layer with rectangular symmetry and the 2R 
crystal consist of two staggered rhombic layers. 

For H2A < H < i?4A, there is a much more com- 
plex cascade of close-packed structures. In the transition 
regime nA — > (n + !)□ for n = 2, 3, on the one hand, 
we recover all of the phases found previously. Here, we 
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FIG. 2. Packing fraction <j> versus dimensionless height H/a. 
The best-packing phases are indicated by symbols on the top 
axis of the middle panel and their packing fractions are shown 
as the full lines. For clarity, the regions below and above the 
best packing are colored differently. The new prism phases 
are denoted by green lines. Dashed and dotted lines denote 
the non-close packed nD [fcc(100)], nA [fcc(lll)], nhcpl, uPa 
and nPn. The top and bottom panel show enlargements of 
the regions where the new prism phases were found. Side 
views (middle panel) and top views (top and bottom panels) 
show the structure of these phases, where white lines denote 
bonds between touching particles. 
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obtain the sequence 2A — > 2Pa — > 2hcpl — ^ 2hcp_L — ^ 
2hcpl -> 2P n -> 3D and 3A -> 3P A -> 3hcpl -> 2P n 
3D, where the following phases are encountered: The 
2n-layered phases uPa and roPn consist of alternating 
prism-like dense-packed n-layered arrays of spheres with 
triangular (A) and square (□) basis shape s f|22l .l23l.l27|). 
Moreover, the 2 n-layered phases rehcp-L 13311 and nhcp- 
like with rectangular symmetry (see [24L l25|) are found. 
For n = 3, however, the nhcp-like phase is only close- 
packed in a tiny regime, whereas nhcp_L is not found at 
all (see Fig.©. 

In the range 3D — > 3 A, the new adaptive prism phase 
IP A is predicted to be close-packed. Representative 
intra-layer touching bonds are indicated by white lines 
in Fig. [2] (upper and lower panel) to underline the sym- 
metry of the corresponding prismatic structure. The IP a 
phase adapts its internal structure flexibly to an increase 
of the slit width H. In fact, the symmetry of its prism 
blocks is rhombic which spans the whole range between 
the square symmetry of the underlying phase 3D and the 
triangular base (see white lines in Fig. [2j upper panel) . 
Likewise, we noticed the stability of the 3Pp prism phase 
with square basis shape (cf. Fig. [2 lower panel) in the 
transition regime 4D — > 4A whose prisms exhibit a longi- 
tudinal shift (ie. in the lenghtwise direction of the prisms) 
in addition to the usual shift perpendicular to the length- 
wise direction of the prisms. Finally, the other densest 
packed phases are multilayered rhombic phases 3P and 
4R as well as a square prism phase 4Pq see Fig. [2j 

To verify our theoretical results, we performed real- 
space experiments with nanometer-sized colloids. We 
employed Polystyrene particles with diameters a in the 
range from 245 nm to 800 nm (Ikerlat Polymers) to study 
a certain of H/a values. We created a confining wedge 
cell with a very small opening angle (10 -4 rad) and slit 
height H < 6um (see |29(). The varying slight height in- 
herent to the wedge geometry allows many transitions be- 
tween different crystals in the same cell. Finally, after the 
sample was dried, we detached the Polystyrene covering 
plate. Some particles stuck to the covering plate during 
its removal resulting in holes in the top layer of particles, 
which allowed us to study the structure in the layers be- 
low. We recorded Scanning Electron Microscope (SEM) 
images from the top facets and side edges by cleaving 
the samples or by Focused Ion Beam milling following 
the crystal planes. 

Concentrating on the regime 3D —¥ 3 A, we found ev- 
idence of the adaptive 2P^ phase. Also, the 3P^ phase 
has been observed for larger plate separation distances. 
As an example, SEM images of 2P A and 3P^ are shown 
in Fig. [3] along with a simulation snapshot for compari- 
son. 

Experiments on colloidal systems, such as ours, are 
necessarily performed at finite pressure. In order to in- 
vestigate the stability of the new prism phases away from 
close packing, we performed Monte Carlo simulations at 



a fixed lateral pressure p = —H~ 1 dF/dA, where F is 
the free energy and A denotes the area of the system. 
This definition of pressure is such that it approaches the 
bulk pressure as H increases. The discovery of new crys- 
tal phases in this and previous theoretical works at infi- 
nite pressure after the previous simulation work that ad- 
dressed the stability at finite pressure begs the question 
how stable these phases are at a high, but finite pres- 
sure We simulated the system at a high pressure 
pa 3 /ksT = 40, for which the system would equilibrate 
within a reasonable time (for comparison the bulk crys- 
tallization pressure is Pa 3, /ksT = 11.56 [HI). The suc- 
cess of cell theory-cffcctivcly a single-particle theory-at 
high densities indicates that phase behavior at high pres- 
sures can be accurately modeled using relatively small 
systems. Our variable-shape simulation box contained 
in x in x n particles, where n is the number of layers 
and 4 < in < 8 [29(. In Fig. [4j we compare finite- 
pressure simulation data to theoretical results at infinite 
pressure. We clearly see that the packing fractions in 
both cases feature a qualitatively similar course. How- 
ever, some phases vanish for finite pressure as this regime 
is dominated by broadened stability regimes of nA and 
nO phases. In detail, the 2P, 3P, 2hcp_L, 3hcpl, 2Pa 
and 3Pa phases are not found for the finite pressure and 
accuracy H/a ± 0.025 chosen in the simulations. As can 
be further seen, the adaptive prism phase 2P4 and 3Pq 
found in this work are stable at this pressure and, there- 
fore, also at all higher pressures. 




FIG. 3. SEM micrographs of the prism phase found in this 
work: 2Pa (a) and 3Pq (c,d). A few particles were removed 
from the top layer upon detachment of the covering plate al- 
lowing access to the structure in the layer below. A simulation 
snapshot, where a particle was also removed (after the simu- 
lation), is shown in (b). White lines indicate the symmetry of 
each phase (a,c) as well as the structure of the prism arrays 
in the side view of 3Pq (d). 
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FIG. 4. Isobaric packing fractions 0, as measured in simu- 
lations at fixed lateral pressure Pia 3 /UbT = 40 for confined 
hard spheres, versus dimensionless separation H/a (empty 
symbols) compared to the theoretical results of Fig. [2] (lines) . 
Each type of empty symbol denotes a different phase. The 
observed phases are indicated by the symbols on the horizon- 
tal lines at the top (theory) and bottom (simulation) of the 
graph. 



close-packed building blocks. 

The adaptive prism phases found here offer new op- 
portunities for several applications. For example, the 
reported structures possess pronounced symmetry direc- 
tions whose alignment can be internally controlled by 
the slit height instead of using external fields (eg. elec- 
tric fields, cf. [36|). As a consequence, these phases can 
serve as switchablc materials. Furthermore, we expect 
an unusual and anisotropic dynamical resp onse of the 
multilaycred prism phases upon shear [37| with possi- 



bly molten grain boundaries which can be exploited to 
tune the rhcological properties of thin crystalline sheets. 
Finally, by varying H , it is possible to tune the whole 
complex cascade of close-packed structures. This may be 
of importance to fabricate nano-sieves or porous mem- 
branes [28[ in a controlled way. 
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PROMETEO/2010/043. F. R.-M. acknowledges the sup- 
port from the EU Marie Curie project APPCOPTOR- 
275150 (FP7-PEOPLE-2010-IEF). 



These simulations help explain the absence of the tri- 
angular prism phase in the experiments (see 27J). We 
also performed simulations with the triangular prism 
phases as initial configuration. At the values for H where 
the triangular prism phase has the highest density of all 
possible phases, the nP& phase appears to consist of n 
only slightly distorted hexagonal layers. At finite pres- 
sure, the small distortions can quickly disappear and a 
regular triangular crystal can be formed. This is a typical 
scenario for crystal-crystal transitions for hard particles, 
where the close packed crystal phase transforms into a 
higher-symmetry crystal with a slightly lower density, 
but a greater entropy, as the pressure is decreased suf- 
ficiently. At larger values of H than those investigated 
here, the triangular prisms are significantly enhanced and 
do not transform so easily to triangular crystals (cf. [23| ) ■ 
Preliminary simulations show that indeed stable triangu- 
lar prim phases can be found for larger values for H (see 
0). 

In conclusion, we explored the close-packed structures 
of hard spheres confined between hard plates in a broad 
range of plate separations by combining theory, exper- 
iment and simulation. We identified adaptive prism 
phases with rhombic symmetry which pack densest in 
certain ranges of the slit width. An adaptive prism phase 
optimizes packing by adjusting its base symmetry flexi- 
bly to the slit width. Also, we showed a high persistence 
of these adaptive prism phases at finite, but large pres- 
sure using experiments and simulations. We anticipate 
that adaptive prism phase will play a key role for even 
higher plate distances, H/a > 3.5, as ideal interpolating 
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